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ABSTRACT 

Recent developments are reviewed in the scaling theory of phase-coherent conduction 
through a disordered wire. The Dorokhov-Mello-Pereyra-Kumar equation for the dis- 
tribution of transmission eigenvalues has been solved exactly, in the absence of time- 
reversal symmetry. Comparison with the previous prediction of random-matrix theory 
shows that this prediction was highly accurate — but not exact: The repulsion of the 
smallest eigenvalues was overestimated by a factor of two. This factor of two resolves 
several disturbing discrepancies between random-matrix theory and microscopic calcula- 
tions, notably in the magnitude of the universal conductance fluctuations in the metallic 
regime, and in the width of the log-normal conductance distribution in the insulating 
regime. 

1. Introduction 

In 1980, Anderson, Thouless, Abrahams, and FisherS proposed a "new method for a 
scaling theory of localization" , based on Landaiier's interpretation of electrical con- 
duction as quantum mechanical transmission.a They considered a one-dimensional 
(ID) chain with weak scattering (mean free path I much greater than the Fermi 
wave length Af), and computed how the transmission probability T scales with 
the chain length L. For L > I an exponential decay was obtained, demonstrating 
localization. In the following decade the scaling theory of ID localization was devel- 
oped in great detail, 3-7 and the complete distribution P(T,L) of the transmission 
probability was found (and hence of the conductance G = T x 2e 2 /h). One can 
thus regard the problem of ID localization as solved, at least for the case of weak 
scattering. 

A real metal wire is not one-dimensional. Typically, the width W is much 
greater than Af , so that the number N of transverse modes at the Fermi level is 
much greater than one. Instead of a single transmission probability T, one now 
has transmission eigenvalues Ti, T2, . . . Tjy. (The numbers T n € [0,1] are the 
eigenvalues of the matrix product ttf , where t is the N x N transmission matrix of 
the wire.) To obtain the distribution of the conductance 
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one now needs the joint probability distribution P(Ti,T2, . . . Tjy,L). This distri- 
bution differs essentially from the distribution in the ID chain, because of strong 
correlations between the transmission eigenvalues. These correlations originate from 
a "repulsion" of nearby eigenvalues. As a consequence of this eigenvalue repulsion, 
the localization length is increased by a factor of N in comparison to the ID case. 
One can therefore distinguish a metallic and an insulating regime. On length scales 
I < L < Nl the conductance decreases linearly rather than exponentially with L. 
This is the (diffusive) metallic regime, where mesoscopic effects as weak localiza- 
tion and universal conductance fluctuations (UCF) occur. The insulating regime of 
exponentially small conductance is entered for wire lengths L > Nl. 

A scaling theory of localization in multi-mode wires was pioneered by Dorokhov,! 
and independently by Mello, Pereyra, and Kumar .i The DMPK scaling equation, 

J = niAj-A,!", (3) 

i<j 

describes the evolution of the distribution function P(Ai, A2, . . . \n, L) in an ensem- 
ble of disordered wires of increasing length. The variables A„ 6 [0, 00) are simply 
related to the transmission eigenvalues by A„ = (1 — T n )/T n . The ensemble is 
characterized by a mean free path I and by a symmetry index (3, which takes on 
the values f , 2, and 4 depending on the presence or absence of time-reversal and 
spin-rotation symmetry (/3 = 2 in the presence of a magnetic field; otherwise, (3 = 1 
or 4 in the absence or presence of spin-orbit scattering). 

The DMPK equation has the form of a diffusion equation in eigenvalue space. 
The function J which couples the degrees of freedom is the Jacobian from the space 
of scattering matrices to the space of transmission eigenvalues. The similarity to 
diffusion in real space has been given further substance by the demonstration!!^ that 
Eq. (||) holds on length scales 3> I regardless of the microscopic scattering properties 
of the conductor (one-parameter scaling). Ecl (|^) was derived by Dorokhov ,@ (for 
[3 = 2) and by Mello, Pereyra, and Kumar (for [3 = 1, with generalizations to 
[3 = 2, 4 in Refs. 11,12) by computing the incremental change of the transmission 
eigenvalues upon attachment of a thin slice to the wire. It is assumed that the 
conductor is weakly disordered (I 3> Af), so that the scattering in the thin slice can 
be treated by perturbation theory. A key simplification is the isotropy assumption 
that the flux incident in one scattering channel is, on average, equally distributed 
among all outgoing channels. This assumption restricts the applicability of the 
DMPK equation to a wire geometry (L ^> W) , since it ignores the finite time scale 
for transverse diffusion. 

The DMPK equation has been studied extensively for more than ten years. 
The strong coupling of the scattering channels by the Jacobian prevented an exact 
solution by standard methods. The problem simplifies drastically deep in the local- 
ized regime (L 3> Nl), when the scattering channels become effectively decoupled. 
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PichardEj has computed from Eq. (^J) the log- normal distribution of the conductance 
in this regime, and has found an excellent agreem ent w ith numerical simulations. 
In the metallic regime (L -C Nl), Mello and StoncElrEJ were able to compute the 
first two moments of the conductance, in precise agreement with the diagrammatic 
perturbation theory of weak localization and UCF. Moregeneral calculations of the 
weak localization effcctlij and of universal fluctuationaij [for arbitrary transport 
properties of the form A = J2 n /Pn)] have been developed, based on linearization 
of Eq. @) in the fluctuations of the A's around their mean positions (valid in the 
large- N metallic regime, when the fluctuations are small) . None of these approaches 
was capable of finding the full distribution function. The purpose of this paper is 
to review some recent work by B. Rejaei and the authorJiZI in which the DMPK 
scaling equation has been solved exactly for the case j3 = 2. 

2. Random-matrix theory and the 1/8 — 2/15 puzzle 

There existed a special and urgent reason for wanting the full distribution function of 
the transmission eigenvalues. We are referring to a disturbing discrepancy^ between 
the random-matrix theory of UCF and the established diagrammatic perturbation 
theory. In order to appreciate the significance of the recent developments, it seems 
worthwhile to discuss this issue in some "historical" perspective. 

In the sixties, Wigner, Dyson, Mehta, and others developed random-matrix 
theory (RMT) into anowerful tool to study the statistics of energy levels measured 
in nuclear reactions.cij It was shown that the fluctuations in the energy level density 
are governed by level repulsion. Mathematically, level repulsion originates from 
the Jacobian J = Yii<j ~ °f the transformation from matrix space to 
eigenvalue space, which depends on the symmetry of the Hamiltonian ensemble — 
but is independent of the mean level density.Ej Thisuniversality is at the origin of 
the remarkable success of RMT in nuclear physics. Ell The universality of the level 
fluctuations is expressed by the celebrated Dyson-Mehta formulate for the variance 
of a linear statistic A = J2 n a i^n) on the energy levels E n . (The quantity A is called 
a linear statistic because products of different E^s do not appear, but the function 
a(E) may well depend non-linearly on E.) The Dyson-Mehta formula reads 

1 1 f°° 

Va,rA= \ dk\a{k)\ 2 k, (4) 

P n Jo 

where a{k) = J^° oo dEe lkE a(E) is the Fourier transform of a(E). Eq. (^) shows that: 
1. The variance is independent of microscopic parameters; 2. The variance has a 
universal 1 / /3-dependence o n th e symmetry index. 

In a seminal 1986-paper,tia Imry proposed to apply RMT to the phenomenon 
of universal conductance fluctuations, which was discovered using diagrammatic 
perturbation theory by Al'tshulernJ and Lee and StonerJ UCF is the occurrence of 
sample-to-sample fluctuations in the conductance which are of order e 2 /h at zero 
temperature, independent of the size of the sample or the degree of disorder - 
as long as the conductor remains in the diffusive metallic regime. The relationship 
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between the statistics of energy levels measured in nuclear reactions on the one hand, 
and the statistics of conductance fluctuations measured in transport experiments on 
the other hand, was used by Muttalib, Pichard, and StonccJ to develop a random- 
matrix theory of quantum transport. (For a review, see Ref. 27.) The RMT of 
quantum transport differs from the RMT of level statistics in two essential ways. 

(i) The first is that the transmission eigenvalues T n are not the eigenvalues of the 
scattering matrix. Instead they are the eigenvalues of tt' , where the transmis- 
sion matrix t is an N x N submatrix of the 2N x 2N scattering matrix of the 
conductor. It turns out that the repulsion of the variables A n = (1 — T n )/T n 
takes the same form as the repulsion of the energy levels E n . More precisely, 
the Jacobian (||) in terms of the A's has the same form as for level statistics. 
Random-matrix theory is based on the fundamental assumption that all cor- 
relations between the eigenvalues are due to the Jacobian. If all correlations 
are due to the Jacobian, then the probability distribution P(Ai, A2, . . . Ajv) of 
the A's should have the form P oc JYliP(Xi), or equivalently, 

P({A„}) = Ccxp[-/?(5>(A l ,A,)+]Tl/(A0)], (5) 

i<j i 

u(Xi,Xj) = -ln|Aj-Aj|, (6) 

with V — — f)~^\np and C a normalization constant. Eq. (|^) has the form 
of a Gibbs distribution at temperature /3 _1 for a fictitious system of classical 
particles on a line in an external potential V, with a logarithmically repulsive 
interaction u. All microscopic parameters are contained in the single function 
V(A). The logarithmic repulsion is independent of microscopic parameters, 
because of its geometric origin. 

(ii) The second difference is that the correlation function of the A's is not trans- 
lationally invariant, due to the positivity constraint on A. This constraint 
A > follows directly from unitarity of the scattering matrix. In contrast, the 
correlation function in the RMT of energy levels is translationally invariant 
over the energy range of interest. Because of this complication, it could not be 
shown that the universality of the fluctuations is generic for arbitrary linear 
statistics on the transmission eigenvalues. In particular, no formula with the 
generality of the Dyson-Mehta formula (|4|) could be derived. The lack of such a 
general theory was being felt especially since mesoscopic fluctuations in trans- 
port properties other than the conductance (both in conductors and super- 
conductors) became of interest. Examples are the critical-current fluctuations 
in Joseph son junctions)^ conductance fluctuations at normal-superconductor 
interfaces,!^ and fluctuations in the shot-noise power of metals Ei This obstacle 
towards the establishment of universality in the RMT of quantum transport 
was finally overcome in 1993,E3 when a technique was developed to compute 
correlation functions by a method of functional derivatives, which does not re- 
quire translational invariance. The analogue could be obtained of the Dyson- 
Mehta formula for the variance of a linear statistic A = J2 n f(T n ) 011 the 
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transmission eigenvalues: 

VarA=-— / dfc|F(/c)| 2 fctanh(7r/c). (7) 
P n Jo 

The function F(k) is defined in terms of the function /(T) by the transform 

fW = /_> eih/ (lT?)' (8) 

The formula (0) demonstrates that the universality which was the hallmark 
of UCF is generic for a whole class of transport properties, viz. those which 
are linear statistics on the transmission eigenvalues. 

The probability distribution (^|) was justified by a maximum-entropy principle 
for quasi-lD conductors.ESO Quasi-ID means L ^> W Ap. In this limit one can 
assume that the distribution of scattering matrices is only a function of the trans- 
mission eigenvalues (isotropy assumption) . The distribution (^|) then maximizes the 
information entropy subject to the constraint of a given density of eigenvalues. The 
function V(X) is determined by this constraint and is not specified by RMT. 

It was initially believed that Eq. (Q) would provide an exact description in the 
quasi-lD limit, if only V(X) were suitably chosen. Ell However, the generalized Dyson- 
Mehta formula (0) demonstrates that RMT is not exact, not even in the quasi- 
lD limit. If one computes from Eq. (Q) the variance of the conductance (Q) [by 
substituting f(T) = G T, with G = 2e 2 /h], one finds 

VarG/Go^i/T 1 , (9) 

independent of the form of V(X). The diagrammatic perturbation theoryBE! of 
UCF gives instead 

VarG/G = 4/ ? " 1 ( 10 ) 
15 

for a quasi- ID conductor. The difference between the coefficients | and is tiny, 
but it has the fundamental implication that the interaction between the A's is not 
precisely logarithmic, or in other words, that there exist correlations between the 
transmission eigenvalues over and above those induced by the Jacobian. 

The | — discrepancy raised the question what the true eigenvalue interac- 
tion would be in quasi-lD conductors. Is there perhaps a cutoff for large separation 
of the A's? Or is the true interaction a many-body interaction, which can not be 
reduced to the sum of pairwise interactions? This transport problem has a counter- 
part in a closed system. The RMT of the statistics of the eigenvalues of a random 
Hamiltonian yields a probability distribution of the form (^L with a logarithmic re- 
pulsion between the energy levelsr!] It was shown by EfetovU and by Al'tshuler and 
ShklovskiiO that the logarithmic repulsion in a small disordered particle (length L, 
diffusion constant D) holds for energy separations small compared to the Thou- 
less energy E c = TiD/L 2 . For larger separations the interaction potential decays 



6 C. W. J. Beenakker 



algebraically. El As we shall discuss, the way in which the RMT of quantum transport 
breaks down is quite different. 

3. Non-logarithmic eigenvalue repulsion 

The method of solution of the DMPK equation is a mapping onto a model of 
non-interacting fermions, inspired by Sutherland's mapping of a different diffusion 
equation.E.3 The case ft = 2 is special, because for other values of f3 the mapping 
introduces interactions between the fermions. The free-fermion problem has the 
character of a onc-dimcnsional scattering problem in imaginary time, which can be 
solved exactly without great difficulties. The reader who is interested in "how it is 
done" is referred to Ref. 17. In this Brief Review we limit ourselves to presenting 
the solution and discussing its implications. 

The DMPK equation (|J) [with = 2] can be solved for arbitrary initial condi- 
tions. We consider the ballistic initial condition limr^o P = JT^ <S(Aj), appropriate 
for the case of ideal contacts. The solution is given by the square root of the Jaco- 
bian (|^) times the determinant of an ./V-dimensional matrix M. The determinant is 
the Slater determinant of the free-fermion problem. The square root of the Jacobian 
comes from the mapping of the DMPK equation onto the Schrodinger equation. The 
solution is 

P({A n },£) = C(L) J 1/2 |DetM|, (11) 
M nm = J dk cxp(-ifc 2 J L/M)tanh(i7rfc)fc 2m - 1 Pi (ife _ 1) (l + 2A„), (12) 

where C(L) is a A-independent normalization factor. Using an integral representa- 
tion for the Legendre functions P„, the matrix elements ( [l2| ) can be rewritten in 
terms of Hermite polynomials Ebm-i, 

/•OO 

M„ m = c / du exp(-\u 2 Nl/ L)(coshu-l-2\ n y 1/2 R 2m -i(^u v / Nl/ L), 

J arccosh(l+2An.) 

(13) 

where c is another constant which can be absorbed in C{L). 

For N = 1, the Jacobian J = 1 and Det M — Mu, so that Eq. ( |TT| ) reduces to 

/>oo 

P(\,L) = C{L) duexp(-± u 2 Z/L)(coshu-l-2A)- 1/2 u. (14) 

Jarccosh(l+2A) 

Normalization gives C(L) = (2vr)- 1 /2^/ jL ) 3 / 2 exp(-iL/0. This is Abrikosov's 
solutionP of the scaling equation for a ID chain. This solution is /3-independent 

a This solution (|l4|)-qf the ID scaling equation was actually obtained as early as 1959 by Gertsen- 
shtein and VasiFcvJ££|in a paper entitled "Waveguides with random inhomogeneities and Brownian 
motion in the Lobachevsky plane. " This remarkable paper on the exponential decay of radio-waves 
due to weak disorder contains many of the results which were rederived in the eighties fpn the 
problem of ID localization of electrons. 3-7 The paper was noticed in the optical literature, E3 but 
apparently not among solid-state physicists. 
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((3 drops out of Eq. (^) for N = 1). Equation ( pi] ) generalizes the ID-chain solution 
to arbitrary N, for the case (3 = 2. 

The Slater determinant can be evaluated in closed form in the metallic regime 
L <C Nl and in the insulating regime L ^> Nl. In both regimes the probability dis- 
tribution takes the form (|^) of a Gibbs distribution with a parameter-independent 
two-body interaction u(Aj, Aj), as predicted by RMT. However, the interaction dif- 
fers from the logarithmic repulsion of RMT. Instead, it is given byEZI 

u(Aj, Xj) = —\ In |Aj — Ai| — | In |arcsinh 2 A^ 2 — arcsinh 2 Ay 2 |. (15) 

The eigenvalue interaction (p5|) is different for weakly and for strongly transmitting 
scattering channels: u — > — In \Xj — Xi\ for Xi, Xj <C 1, but u — ► — ^ In |Aj — \ \ for 
Xi,Xj 3> 1. For weakly transmitting channels it is twice as small as predicted by 
considerations based solely on the Jacobian, which turn out to apply only to the 
strongly transmitting channels. 

In the metallic regime L <C Nl, the method of functional derivatives of Ref. 18 
can still be used to compute the variance of a linear statistic, since this method 
works for any two-body interaction. Instead of Eq. (ffl), one now obtains for the 
variance the formula 

^A=i\r dk _ (16) 

P n Jo 1 + cotanh(i7rfc) 

F(k)= Tdx^ff^-). (17) 
J — oo \ cosn xj 

This result was obtained for f3 = 2 from the exact solution given above,0 and 
independently for all (3 £ {1,2,4} by the perturbative method of Chalker and 
Macedo.li3 Substitution of f(T) = T now yields ^ instead of | for the coefficient of 
the UCF, thus resolving the discrepancy between Eqs. (||) and (|l0|). The conclusion 
is that the discrepancy with RMT originated from a reduced repulsion of weakly 
transmitting channels. 

In the insulating regime L ^> Nl, all A's are exponentially large, and the inter- 
action ( p"5| ) may be effectively simplified by u(Xi,Xj) = — |ln|Aj — Aj|. This is a 
factor of two smaller than the interaction (^) predicted by RMT. This explains the 
factor-of-two discrepancy between the results of RMT and of numerical simulations 
for the width of the log-normal distribution of the conductance £3 RMT predicts 
Var lnG/Go = —(IhG/Gq), which is twice as small as the result 

Var In G/G = -2 (In G/G ) (18) 

which follows from the exact solution of the DMPK equation for /3 = 2. As shown by 
Pichardjll the relationship (|l8| ) between mean and variance of In G/Gq remains valid 
for other values of f3, since both the mean and the variance have a 1//3 dependence 
on the symmetry index. 
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4. Outlook 

We conclude by mentioning some directions for future research. So far only the case 
j3 = 2 of broken time-reversal symmetry has been solved exactlyO In that case the 
DMPK equation (^) can be mapped onto a free-fermion problem. For /? = 1,4 the 
Sutherland-type mapping is onto an interacting Schrodinger equation. It might be 
possible to solve this equation exactly too, using techniques developed recently for 
the Sutherland Hamiltonian.Ell From the work of Chalker and Macedoli3 we know 
that the two-point correlation function of the eigenvalues in the large- N limit has a 
simple 1 / /3-dependence on the symmetry index. This poses strong restrictions on a 
possible /3-dependence of the eigenvalue interaction, which can only differ from the 
form ( |l5|) derived for (3 = 2 on intervals A A ~ L/Nl <C 1 comparable to the spacing 
between the A's. 

It might be possible to come up with another maximum-entropy principle, differ- 
ent from that of Muttalib, Pichard, and Stone,cl which yields the correct eigenvalue 
interaction (|l^) instead of the logarithmic interaction (|^) . Slevin and NagacH have 
proposed an alternative maximum-entropy principle, but their distribution function 
does not improve the agreement with Eq. (|l^). It would be particularly worthwhile 
to find an intuitive explanation for the halving of the logarithmic interaction for 
weakly transmitted scattering channels. 

To go beyond quasi-one-dimensional geometries (long and narrow wires) remains 
an outstanding problem. A numerical study of Slevin, Pichard, and MuttalibH has 
indicated a significant break-down of the logarithmic repulsion for two- and three- 
dimensional geometries (squares and cubes) . A generalization of the DMPK ■yma- 
tion (||) to higher dimensions has been the subject of some recent investigations e3e3 
It remains to be seen whether the method reviewed here for Eq. (Q) is of use for 
that problem. 
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